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We apply majorization theory to study the quantum algorithms known so far and find that 
there is a majorization principle underlying the way they operate. Grover's algorithm is a neat 
instance of this principle where majorization works step by step until the optimal target state 
is found. Extensions of this situation are also found in algorithms based in quantum adiabatic 
evolution and the family of quantum phase-estimation algorithms, including Shor's algorithm. 
We state that in quantum algorithms the time arrow is a majorization arrow. 
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Majorization is the natural ordering on probability dis- 
tributions. One probability distribution is more uneven 
than another one when the former majorizes the latter. 
Furthermore, majorization implies an entropy decrease, 
thus the ordering concept introduced by majorization is 
more restrictive and powerful than the one associated to 
the Shanon's entropy. The goal of this work is to show 
that all known efficient quantum algorithms obey a ma- 
jorization principle, in a way to be made precise later. 

The classical theory of majorization was first intro- 
duced by Muirhead Q and later developed by Hardy, 
Littlewood and Polya in their study of symmetric means 
Majorization was early studied by economists in the 
beginning of the twentieth century in order to formal- 
ize the concept of unevenness in the distribution of in- 
come. In 1905, Lorenz pointed out that one distribu- 
tion can be said to be more uneven than another pre- 
cisely when it majorizes the other ||. Likewise, Dalton 
in 1920 stated his principle of transfers showing that a 
distribution is less uneven than another if it can be ob- 
tained from the other by transferring some income from a 
richer to a poorer income-receiver. Moreover, majoriza- 
tion has found many applications in classical computer 
science like stochastic scheduling, optimal Huffman cod- 
ing, greedy algorithms, etc. 

In quantum information theory, majorization charac- 
terizes when two quantum bipartite pure states can be 
connected via Local Operations and Classical Communi- 
cation ^, D . This result shows that this connection is in- 
deed possible when there exists majorization between the 
vectors of eigenvalues (weights) of the partial von Neu- 
mann entropies associated to each bipartite state. A fur- 
ther application of majorization in quantum information 
theory corresponds to the problem of Hamiltonian simu- 
lation There, strong restrictions based on majoriza- 
tion theory limit the possibility to simulate a proposed 
quantum evolution from a different given Hamiltonian 
complemented with local unitary transformations. Ma- 



jorization is also present in quantum measurement theory 
and in the separability problem. 

Majorization is often defined as a binary relation de- 
noted by -i on vectors in R d . We need to fix notations 
by introducing some basic definitions. 
Definition 1. For x, y G K d , 

x^yiffjPr^^r^' fc = 1 >-< d - 1 (!) 
[2^=i ^[i] = L<=iV[i]. 

where [z^j . . . z^]] := sort^(z) denotes the descendingly- 
sorted (non-increasing) ordering of z G M. d . An imme- 
diate consequence is that majorization is a partial order 
for sorted vectors in R d . 

Definition 2. If it exists, the least element x\ (greatest 
element x s ) of a partial order like majorization is defined 
by the condition x\ -< x,\/x G R d (x -< x g ,\/x G R d ). 

In this letter we address the following basic prob- 
lem of elucidating what is the role, if any, played by 
majorization in the way quantum algorithms operate. 
We find, indeed, that there is a majorization principle 
underlying the way quantum algorithms work that we 
shall now state more precisely. Let us denote by l^m) 
the pure state representing the state of the register in 
a quantum computer at an operating stage labeled by 
m = 0, 1, . . . , M — 1, where M is the total number of 
steps of the algorithm. We can associate naturally a set 
of sorted probabilities [py], x = 0, 1, . . . , 2" — 1 to this 
quantum state of n qubits in the following way: decom- 
pose the register state in the computational basis i.e., 

I*m) := Ex=o lc *K> with i\ x ) '■= koici • • • aJn-i}}^ 1 
denoting the basis states in digital or binary notation, 

respectively, and x := Y^ijZo x j^- The sorted vec- 
tors to which majorization theory applies are precisely 
\P\x\] [l c [:r]| 2 ]- Thus, in quantum algorithms we shall 
be dealing with probability densities defined in with 
d = 2 n . With these ingredients, our main result can be 
stated as follows: in the quantum algorithms known so 
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far, the set of sorted probabilities [p^] associated to the 
quantum register at each step m are majorized by the 
corresponding probabilities of the next step 



m+li 




(2) 



This is a strong result for it means that majorization 
works locally in quantum algorithms. i.e., step by step, 
and not just globally (for the initial and final states). Our 
starting point is the majorization analysis of Grover's 
algorithm @. 

Grover's algorithm. This quantum algorithm solves ef- 
ficiently the problem of finding a target item in a large 
database. The algorithm is based on a kernel that acts 
symmetrically on the subspace orthogonal to the solu- 
tion. This is clear from its construction 



K := U.TJ, 
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2\s) 



1, Uy. 
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where \s) := \ x ) an< ^ \Vo) is the searched item. 

Theorem: The set of probabilities to obtain any of the 
N possible states in a database is majorized step by step 
along the evolution of Grover's algorithm when starting 
from a symmetric state until the maximum probability 
of success is reached. 

Proof. To prove this result we write [pu] as the set of 
sorted probabilities of finding the state \x) when per- 
forming a measurement. We call [pLi] the set of sorted 
probabilities after one single application of Grover's ker- 
nel. The theorem is equivalent to prove that [p[ x ]] -< [pL]] 
until pi, the probability of finding the correct solution, 
reaches its maximum value. 

The hypothesis of symmetry imposes that the proba- 
bilities of finding each of the N outputs at some point 
during the implementation of Grover's algorithm can be 
ordered in the list 



i-p i-p 

r> N-l ' N-l ■ 



l-p 

N-l 



(4) 



where p is the one associated to the correct output. After 
one further action of the kernel these probabilities will be 



f ' N-l ' N-l ' 



i-v' 

N-l 



(5) 



We first need to prove that Grover's algorithm increases 
the probability of success monotonically, that is p 1 > p, 
till it reaches a maximum and then decreases also mono- 
tonically This part of the proof relies on the fact that 
the Grover algorithm can be described in a reduced two- 
dimensional space f|],[[)), which follows from the sym- 
metry of the subspace orthogonal to |yo)- In this case, 
the dynamics can be reduced to a two-state system, 



{\yo), \y_o }}• Grover's kernel on this space acts as a rota- 
tion lied 
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( cos 8 — sin 9 
ysind cos 9 



(6) 



where cos ( 



1 — j£ . Starting from the symmetric state 



Jn y 1 n ) ' 



m applications of the kernel lead to 



(7) 



4= cos m6 — \ / 1 — 4? sin md \ 

I ^ V ^ . (8) 

-4= sin mO + J 1 — cos md J 

The projection onto the upper component corresponds 
to the probability amplitude which, thus, evolves mono- 
tonically until it reaches a maximum. 

Returning to the original problem, we can now check 
that all probabilities evolve in such a way that majoriza- 
tion works smoothly: 

P < P', 

(N - 2)p + 1 (N - 2)p' + 1 



N - 1 



N - 1 



(9) 



(N — m — l)p + m (N — m — \)p' + m 



N- 1 



N- 1 



Thus [p[ x ]} -< \pi,] and @ holds true. ■ 

Majorization works in a simple way in Grover's algo- 
rithm. Nevertheless, the proof does not hold when the 
initial distribution of probabilities is not symmetric in the 
subspace orthogonal to the solution. It is indeed easy to 
find numerical counterexamples to the majorization prin- 
ciple in absence of symmetry. We realize that this corre- 
sponds to starting with a quantum state |s) whose set of 
probabilities is the least element of the majorization we 
have introduced to study quantum algorithms. We shall 
see that this fact also happens in the rest of algorithms 
below. 

Quantum adiabatic evolution algorithms. Grover's algo- 
rithm can be mapped onto the evolution of the homoge- 
neous state |s) into the solution |0) driven by a simple 
Hamiltonian Farhi et &1. have proposed to use the 
adiabatic evolution to guarantee that the system remains 
in the fundamental state and reaches the target solution 
in the end [0. More precisely, the idea consists of set- 
ting up a Hamiltonian of the form 



H (10 — (l ~ TJ Ho + 



(10) 



such that |s) is the ground state of Ho and |0) is the 
ground state of H\ . For large enough T, the evolution will 
be adiabatic and the system will remain in the ground 
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state all along the flow. The adiabatic theorem dictates 
that T must scale as the inverse squared of the mini- 
mum gap of the system. The question we address here is 
whether this evolution respects majorization. 

Although the system contains n qubits, 2" possible 
states, the adiabatic evolution can be computed using a 
subspace if sufficient symmetry is present. The simplest 
example is to consider the Hamiltonian 

ff(±)=-| s )( s |(l-i)-|0)(0|| (11) 

and the initial state |s). In this particular case, the evo- 
lution can be computed using a reduced two state Hilbert 
space. More precisely 

|s) = -^(|0)+V2^I|0 ± )) (12) 

Then the Hamiltonian written in the basis {|0), jO -1 "} 
reads 

g (*)=-(!-*) (vfc 1|f)-?(<S 2) ( i3 ) 

It is possible to verify numerically that when T ~ 4 2™ 
the probability follows the graphic shown in Fig. [p. An 
argument similar to the previous theorem indicates that 
symmetry imposes majorization for the complete set of 
probabilities. Shorter T lead to evolutions that do not 
hit the solution with probability one, while a larger T 
smooths this evolution. Once the maximum is attained, 
the probabilities oscillate and majorization is obviously 
lost. 
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FIG. 1: Evolution of the probability of finding the target 
state (bold) and other states (dashed) for n = 6. 

It is worth mentioning that a combination of Hq and 
Hi chosen as above but mixed with no time dependence 
leads to a Hamiltonian that rotates the ground state in 
the manner of the previous theorem. Then, the solution 
is obtained in T = ^2% with probability 1. This is pre- 
cisely the scaling law found in Grover's algorithm. 

A more refined test for the majorization principle cor- 
responds to the Hamiltonian evolution proposed by Farhi 




FIG. 2: Curves for p\, pi + pi and p\ + P2 + P3 for 
n = 4. The Failure of majorization (monotonicity) for 
fast evolution, T — 4 2™, in the upper curves goes away 
for slower evolution, T = 7 2". 

ct al. as a natural starting point for any adiabatic evo- 
lution [fn]. Let us consider the following choice 

H = £(1 (14) 

i— l,n 

This Hamiltonian acts as an eraser of information and 
has the state \s) as its ground state. Furthermore, it 
allows for a decomposition of the Hilbert space into n + 1 
symmetric subspaces. Finding the target instance |0) 
amounts to solving the dynamical evolution in this (n + 
l)-dimcnsional Hilbert space. Let us denote as \k) as the 
symmetric space with k qubits in the state |1) and the 
rest in |0). The Hamiltonian becomes 

H = ^I-N, (15) 

where the elements of the symmetric matrix N are given 
by 

A numerical solution of the evolution is now easy to per- 
form. For T > 7 2™, the system indeed evolves along the 
ground state and majorization holds for the set of n + 1 
probabilities, as shown in Fig. ||. Shorter evolutions per- 
form poorly and fail to verify the majorization principle. 
We conclude that quantum algorithms based on adia- 
batic evolution naturally fulfill a majorization principle 
provided that the Hamiltonians and initial state are cho- 
sen with sufficient symmetry and the evolution is slow 
enough. 

Quantum phase- estimation algorithms. These represent 
a large family of quantum algorithms that include as par- 
ticular instances the order-finding problem, Shor's algo- 
rithm |Q, discrete logarithms, etc. The basic prob- 
lem is: given an arbitrary unitary operator U and one 
eigenvector \v), estimate the phase <f> of the correspond- 
ing eigenvalue U\v) := e~ 27r '*|w), (f> e [0, 1), with n bits of 
accuracy. The efficient quantum solution of this problem 
can be encoded in the quantum circuit shown in Fig. ||, 
and we shall always refer to this circuit when perform- 
ing the majorization analysis stepwise. The algorithm 
has clearly two parts: i) application of Hadamard gates 
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Ujj and controlled-?/^ gates, j = 0, 1, . . . , n — 1; ii) ap- 
plication of the quantum Fourier transform (QFT) Uf- 
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FIG. 3: a) Quantum circuit implementing the phase- 
estimation algorithm constructed from Hadamard gates 
U H , controlled-?/ gates acting as |0)(0| (g) 1 + |1)(1| ® U, 
and the QFT. Dashed lines represent time steps for ma- 
jorization testing, b) An example of QFT decomposition 
into elementary gates for n = 3 qubits. 

Part i). The whole quantum register is made up of 
first and second registers. The initialization stage is 
such that the quantum computer is in the state l^m) := 
1 00 . . . 0}|v), where the first register has been prepared at 
the state |0) for short, and the second holds the eigen- 
vector of U. In what follows, we denote by [pS] the 
sorted probabilities distributions of the first register, at 
time steps m = 0,l...,n+l that we show in Fig. || as 
time slices. 

Clearly, the probability distribution of |^i n ) is a great- 
est element of the majorization. However, an ap- 
plication of the Hadamard gates yields a lowest ele- 
ment as in Grover's algorithm. Thus, our starting 
point for majorization is |* ) : = (U§ n ® l)|*i n ) = 
S-^E'lo 1 \ x )\ v )- Th cn, \p° x] ] = [2-"],Vx. 

Next, a series of controlled- U 2,1 gates encompassing 
time steps from t% to t n (Fig. ||) arc applied. The out- 
come of these steps is the factorized state 



|*„) =2- n/2 [\Q) + er 2 ™ 2 " ^|1)] 

2 n -l 
=2 -n/2 ^ 



e- 27rix,p \x)\v) 



x=0 



'11)] 



(17) 



As the action of these gates only introduces phases locally 
in the computational states, then we obtain again the 
uniform distributions [pS] = [2 _n ], Va;, m = 0, 1, . . ., n. 



Part ii). Although the local phases in \^f n ) do not play 
any role in majorization, so far, they become relevant 
when combined with the application of the QFT on the 
first register, due to interference of quantum amplitudes. 
The state after time step t n+ i (Fig. ||) is 



|* n +i) := (E7>®l)|tf n >=2- 



2 n — 1 

E e 

X :V— 



-2vix(<j>-y /2 n ) 



Now, p 



n+1 



El 



\y)\v). 

(18) 

-i e -27r t2; (0-y/2")|2 ma j 0I . izes trie 

least element distribution at step m = n. Interestingly 
enough, there is a stronger majorization working step- 
wise when the QFT is applied by means of its canonical 
decomposition in terms of n Hadamard and n(n — l)/2 
controlled-phase gates |14"| . For concreteness, we show 
such decomposition in Fig. |^b) for n = 3 qubits and with 
the corresponding time slices (majorization checkpoints). 
The proof of this result relies on the recursive application 
of the following inequalities 
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(lie 



2-Kia±(y,4>) 




(19) 



where, at each step, a± depends on y, <p in a computable 
way fl5|| . To illustrate this fact, we show in Fig. I| a nu- 
merical plot for n — 3 qubits in the form of a Lorenz 
diagram: partial probability sums vs. x, for each time 
step. Therefore, as a consequence of our analysis we 




FIG. 4: Lorenz diagram (partial probability sums) for 
the quantum phase-estimation algorithm with <p = 0.2 
and n = 3 qubits as in Fig. ||. It shows how majorization 
works along the time arrow Q ^ O (} ^ A. 

find that the majorization principle is working locally 
in algorithms like order-finding a r = 1 mod A, where 
the unitary operator is given by U\x) := \ax mod N) 
and (j) = 1/r; Shor's algorithm, where order-finding is 
used combined with controlled-?/ gates implementing the 
modular exponentiation; Chuang's algorithm for quan- 
tum clock synchronization, where U := U cno tUTQpU cno t 
and Utqp is the so called Ticking Qubit Protocol flT6| ; 
etc. 

Conclusion. Efficient quantum algorithms are scarce 
as compared with their classical counterparts, suggest- 
ing that we are missing the basic principles for quantum 
algorithm design ]l7j ]. In this note, we have produced 



evidence for the general idea that there is a majorization 
principle acting step by step during the time evolution [6 
in efficient quantum algorithms. We may say that ma- [7 
jorization is a sort of driving force for such algorithms. [8 
Learning to tame majorization may be useful for devis- 
ing quantum algorithm design. When majorization is not [9 
at work, the quantum algorithm is neither efficient nor 
successful. [10 
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